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ABSTRACT 

We consider the neutron star (NS) of the magnetar type inside the massive binary 
system. We determine the conditions under which the matter from the stellar wind 
can penetrate the inner magnetosphere of the magnetar. At some distance from the NS 
surface, the magnetic pressure can balance the gravitational pressure of the accreting 
matter creating very turbulent, magnetized transition region. It is suggested that this 
region provides good conditions for acceleration of electrons to relativistic energies. 
These electrons lose energy on the synchrotron process and the Inverse Compton (IC) 
scattering of the radiation from the nearby massive stellar companion, producing high 
energy radiation from the X-rays up to ~TeV 7-rays. The primary 7-rays can be farther 
absorbed in the stellar radiation developing the IC pair cascade. We calculate the 
synchrotron X-ray emission from primary electrons and secondary e ± pairs and the 
IC 7-ray emission from the cascade process. It is shown that the quasi-simultaneous 
observations of the TeV 7-ray binary system LSI +61 303 in the X-ray and the TeV 
7-ray energy ranges can be explained in such an accreting magnetar model. 
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1 INTRODUCTION 

The accretion process can occur onto the neutron star in the 
case of relatively slowly rotating objects even if the neutron 
stars are characterized by a super-strong surface magnetic 
field (so called magnetars). Up to now, magnetars have not 
been observed inside the binary systems. However, very re- 
cently it was claimed that X-ray emission characteristic for 
magnetars has been observed from the object LSI +61 303, 
which belongs to the rare class of massive binary systems 
emitting TeV 7-rays (Dubus & Giebels 2008). The observed 
outburst of thermal X-ray emission might be due to the tran- 
sient accretion of matter onto the surface of the neutron star. 
Moreover, quasi-periodic X-ray emission, with the frequency 
of 0.2 Hz, has been reported from the direction of this bi- 
nary (Ray, Hartman et al. 2008). Whether this emission is 
in fact related to this massive binary, it is at present an open 
issue (see Rea & Torres 2008). However, if it is true, then it 
makes sense to speculate that the compact object inside the 
LSI +61 303 is an accreting magnetar. 

The aim of this paper is to investigate such a hypothe- 
sis. We apply the general scenario for the 7-ray production 
similar to that one recently proposed for accreting classical 
and millisecond neutron stars (Bednarek 2008). It is shown 
that also accreting magnetars can accelerate electrons in the 
magnetized, turbulent region inside the NS magnetosphere 



when the gravitational pressure of the accreting matter is 
balanced by the magnetic pressure of the rotating NS mag- 
netosphere. We calculate the high energy radiation produced 
in the synchrotron and inverse Compton processes by pri- 
mary electrons and secondary e cascade pairs immersed in 
the magnetic fields of the magnetar and the massive star, 
and in the thermal radiation from the massive star surface. 



2 DESCRIPTION OF THE MODEL 

We consider a compact binary system containing rotating 
strongly magnetized neutron star (magnetar) and a massive 
companion of the 0,B type star. It is assumed that a mass 
from the stellar wind is effectively captured by the strong 
gravitational potential of the neutron star (NS). Depending 
on its period and surface magnetic field, the accretion pro- 
cess onto NS can occur in different phases. According to the 
classification scheme of Lipunov (1992), the accretion pro- 
cess can occur in the phase of accretor (for relatively slow 
rotators) or in the phase of propeller. The accretor phase in 
the case of classical and millisecond pulsars has been recently 
discussed by Bednarek (2008). In this paper we concentrate 
on the propeller phase for the strongly magnetized neutron 
stars, i.e. magnetars. In this case, the matter from the stel- 
lar wind can penetrate below the light cylinder radius of the 
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rotating NS magnetosphere. This matter extracts rotational 
energy from the NS as a result of interaction of the free 
falling matter with the rigidly rotating inner NS magneto- 
sphere. As a result of this interaction, a very turbulent and 
magnetized transition region is formed. We consider possible 
acceleration of electrons in such a medium. 

The distance at which the magnetic field starts to dom- 
inate the dynamics of the in-falling matter (the AlfVen ra- 
dius) can be estimated by comparing the magnetic field en- 
ergy density with the kinetic energy density of the wind, 



B 2 A /87v = pv 2 { /2, 



(1) 



where B\ is the magnetic field in the inner neutron star 
magnetosphere, p = M acc /(ivR\v{) is the density of the ac- 
creting matter, V{ — (2GMns/Pa) 1,/2 is the free fall velocity 
of the accreting matter, Pa is the Alfven radius, and G is the 
gravitational constant. The matter in this transition region 
is very turbulent and strongly magnetized providing good 
conditions for acceleration of particles to high energies. Let 
us estimate the location of this region from the surface of 
the neutron star. By applying Eq. [1] and assuming that the 
magnetic field in the neutron star magnetosphere is of the 
dipole type, i.e. Ba = Bns(Pns/Pa) 3 , we obtain 



Ra 



5.5 x 10 9 Bfi 7 M 16 



2/7 



(2) 



where the magnetic field at the neutron star surface is Bns = 
10 14 Bi4 G and the accretion rate is M = 10 16 A/i 6 g s" 1 
(see also Baan & Treves 1973). Then, we can estimate the 
magnetic field strength at the transition region, 



Ba = 6WM°''B 



6/7 R -5/7 



(3) 



The accretion according to the propeller scenario can oc- 
cur provided that some conditions are fulfilled. At first, the 
radius of the transition region has to lay inside the light 
cylinder radius of the neutron star, i.e Ra < Plc ~ cP/2tt, 
where P — lPi s is the rotational period of the neutron star, 
and c is the velocity of light. The above condition is fulfilled 
when, 



Pi > l.2Bt{ 7 M" 2/7 . 



(4) 



At second, the rotational velocity of the magnetosphere 
at Ra has to be longer than the keplerian velocity of the 
accreting matter. The rotational velocity, 

v rot = 2ttP a /P « 3.5 x W 10 Bf /7 ' 



' 1 4 -^16 



2/7 /Pi cms" 1 . (5) 



is larger than the keplerian velocity, 



v k = (GA/ ns /Pa 



,V2 



1.8 x 10 s B- 



-2/7 M l/7 



for the NS with the rotation periods, 



Pi < 192P 



6/7, 



-3/7 



(6) 



(7) 



This last condition on the NS period separates the propeller 
phase (lower periods) from the accretor phase (larger peri- 
ods). 

The third condition relates the Alfven radius, Pa, to 
the capturing radius of the matter from the stellar wind by 
the NS. It is determined by the balance between the kinetic 
energy of wind with its potential energy in the gravitational 
field of the NS, 



Pc 



2GM NS /vl 



3.7 x W 10 vo 



cm, 



(8) 



where « w = 10 8 «8 cm s _1 is the velocity of the surrounding 
matter measured in respect to the NS. The accretion from 
the stellar wind occurs when P c > Pa, which happens for, 



Pi4 < 28«~ 7/2 AP\ /2 . 



(9) 



This condition is not restrictive even for magnetars in the 
case of the accretion rates considered below. 

In the conditions expected for the transition region 
(strongly magnetized and very turbulent medium), parti- 
cles should be efficiently accelerated. In this paper we con- 
sider only acceleration of electrons. The acceleration rate of 
electrons with energy E (and the Lorentz factor 7) can be 
parametrized by, 



: fcP/r L » 970£-iM 1 6 6 /7 B M 5/7 



erg s , (10) 

where £ = 10 -1 £-i is the acceleration parameter, c the ve- 
locity of light, tl = E/eBA the Larmor radius, and e elec- 
tron charge. During the acceleration process electrons suffer 
energy losses mainly on the synchrotron process and on the 
inverse Compton scattering of radiation from the massive 
star. 

Electrons lose energy on IC process in the Thomson 
(T) and the Klein-Nishina (KN) regimes. Let us estimate 
the photon energy densities from the massive star (p*) in 
the location of the NS, 

^ = 4^f^y w61xl0 3 T 4^y ergcm _ 3) 



D 



where T* = 3 x 10 4 T4 K and a is the Stefan-Boltzmann 
constant, and D is the distance of the NS from the massive 
star. 

Let us also estimate the energy density of the magnetic 
field at the transition region (see Eq. [3}, 



PB 



Bi 



1.4 x 10 4 Mi 12/7 B 14 



10/7 -3 
erg cm , 



(12) 



The energy losses for each process (synchrotron and IC in 
the T regime) can be calculated from, 



(4/3)ccttP7 



2.7 x 10- 14 p {cap , B , i)7 2 ergs" 1 , (13) 



where or is the Thomson cross section. 

The energy losses of electrons at large energies are dom- 
inated by the synchrotron process since the IC losses in the 
radiation field of the massive star decline due to the Klein- 
Nishina cross section. Therefore, the maximum energies of 
accelerated electrons are determined by the balance between 
energy gains (Eq. HOf) and energy losses on synchrotron pro- 
cess (Eq. I13p . They can be expressed by, 



In 



1.5 x lO^B^X"/ 77 . 



(14) 



It is clear that for the realistic parameters of the model 
electrons can reach TeV energies. 

The maximum power available for acceleration of elec- 
trons is limited by the energy extracted from the rotating 
neutron star by the in-falling matter. This matter from the 
stellar wind is accelerated to the velocity of the magnetic 
field lines at Pa- The power which has to be transfered from 
the rotating NS to the matter can be estimated from 

Lace = Macc« 2 o t /2 « 6 x 10 36 B^fM^Pf 2 erg s~\ (15) 

By Using Eq. [4] we can estimate the maximum power which 
can be extracted in terms of such a model in the propeller 
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stage, 

Lace « 4 x 10 36 Mi 6 erg s' 1 . 



(16) 



A part of this maximum possible power, rj, can go on accel- 
eration of electrons. 

The total rotational energy of the magnetar is, 



P rot = Iu 2 /2 w 4.4 x 10 46 Pf 



erg. 



(17) 



where 7 is the angular momentum of inertia of rotating NS 
calculated for its mass equal to 1.4M Q and the radius of 10 
km, and ui = 2n/P is its angular velocity. This energy reser- 
voir is enough to provide sufficient power to the accreting 
matter for a time of the order of thousand of years. Note that 
a neutron star in a specific binary system can come through 
the cycles of rotational energy losses (during the propeller 
phase) and energy gains (during the accretor phase). These 
phases depend on the accretion rate from the massive star 
which can change due to the stellar variability or evolution. 
So then, the angular velocity of magnetars may also increase 
at some periods as in the case of millisecond pulsars. 

For the known parameters of the binary system, we can 
estimate the accretion rate onto the magnetar, 

M acc = R 2 C M,/{AD 2 ) « 5 x 10 15 A/_ 7 /(« 8 D 12 ) 2 g u -1 . (18) 

where the capturing radius is given by Eq. [5] D = 10 12 Di2 
cm is the separation of the stars inside binary system, 
M* — 10~ 7 M_7 Mq yr _1 is the mass loss rate of the 
massive star, and w w = 10 s vg denotes the wind velocity 
and/or the NS velocity on its orbit around the massive star 



0.377X 



-1/2 



For the luminous stars, the wind velocity, 



of the order of ~ 10 km s~ , is usually larger than the 
NS velocity. The B type stars, present inside the TeV 7-ray 
binaries, have typical mass loss rate in the range 10 - ' 6-7 - 1 
Mq yr _1 and stellar wind velocities ~ (1 — 3) x 10 3 km 
s _1 . Then, the typical accretion rates onto the NS are in the 
range ~ 10 13 ~ lr g s _1 , for the separation of the companion 
stars below 10 radii of the massive star. 



3 X-RAY AND GAMMA-RAY EMISSION 

In this model, radiation is produced by electrons in the 
synchrotron and IC process. We assume that primary elec- 
trons in the transition region reach the power law spectrum, 
oc E~ s up to -Bmax = ?7i c 7max- The energy losses of elec- 
trons are dominated at the highest energies (close to -E max ) 
by the synchrotron process occurring in the magnetic field of 
the transition region. The 7-ray production at these highest 
energies in the IC scattering of thermal radiation from the 
massive star occurs typically at a lower rate than the syn- 
chrotron process but at lower energies IC process can domi- 
nate. The IC 7-rays have large enough energies that in some 
cases they can initiate IC e pair cascades in the radiation of 
the massive star. Such cascade can develop provided that the 
optical depths for electrons in the anisotropic radiation of 
the massive star are large enough. These optical depths have 
been calculated in the most general case of the anisotropic 
radiation of the massive star by e.g. Bednarek (1997, 2000). 
They can be easily scaled for the case of the massive star 
with the arbitrary parameters (see Bednarek 2008). 

We developed the Monte Carlo code which is able to cal- 
culate the synchrotron radiation from the primary electrons 



(magnetic field of the transition region), the secondary cas- 
cade e pairs (in the magnetic field of the massive star), and 
the 7-ray spectra produced in IC e pair cascade occurring 
in the anisotropic radiation of the massive star. This code 
base on the earlier code discussed for the binary systems (e.g. 
Bednarek 2000). In order to calculate the synchrotron spec- 
trum produced by e ± pairs inside the volume of the binary 
system, we have to define the distribution of the magnetic 
field around the massive star. The strength of the magnetic 
field as a function of distance from the center of the massive 
star can be described by the following equations (Usov & 
Melrose, 1992), 



B(r) ~ x 



R 3 J(R*, A r 2 



R*,a < r, 



(19) 



where B* is the surface magnetic field of the massive star, 
i?*,A is the Alfven radius within the massive star wind which 
can be derived by solving the equation, (1 — -R*/r*,A) = 



TV A) 



It has simple approximate solution, 



R t 



R* x 



1 + 6 

e r/4 ) 



?«1 



(20) 



where £ = B* 2 R- k /(Mv cx> ), M is the mass loss rate, and Doo 
is the stellar wind velocity. This approximation for the mag- 
netic field structure is satisfactory inside the most important 
part of the binary system (i.e. within ~ 10i?*) for the stars 
which typical rotation velocity is lower than ~ 10% of their 
wind velocity. Note that the larger values of the surface mag- 
netic field of the massive star do not necessary mean larger 
values of the magnetic field at some distance from the stellar 
surface (as would be expected from the simple analysis of 
Eq. 1 19 j) since B(r) depends also on the Alfven radius which 
in turn depends on the surface magnetic field and the mass 
loss rate of the massive star. 



4 THE CASE OF LSI +61 303 

As an example, we apply the model for the massive binary 
system LSI +61 303 which has been discovered as a TeV 7- 
ray source by the MAGIC Collaboration (Albert et al. 2006) 
and confirmed by the VERITAS Collaboration (Acciari et 
al. 2008). The TeV 7-ray emission from this binary system 
is observed only at a specific range of phases close to the 
apoastron. No any TeV 7-ray emission (above ~ 200 GeV) 
has been observed up to now close to the periastron passage. 
The multi-wavelength observations of LSI +61 303 show a 
hint of correlation between X-ray and TeV 7-ray emission 
(Albert et al. 2008). LSI +61 303 belongs to the rare class of 
sources observed in the whole high energy range from X-rays 
up to TeV 7-rays (Tavani et al. 1996, van Dijk et al. 1996, 
Kniffen et al. 1997). 

The massive star in this binary system belongs to a class 
B0 V. It is characterized by the surface temperature T* = 
29850 K and the radius R* — 6.7Rq (see e.g. Harmanec 
1988). The basic parameters of the binary system are not 
exactly known. Two sets of parameters have been proposed 
by Casares et al. (2005) and Grundstrom et al. (2007). The 
semimajor axis, a, the eccentricity, e, the angle of periastron 
passage, w p , the phase of the observer, u> b B , the inclination 
angle, a, the distance at periastron and apoastron, r p and 
r a , are shown in Table 1. 
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Table 1. Basic parameters of the models. 



asini u) v e ^obs a r p r a M* ihr 

model I (Grundstrom et al. 2007) 8.4 R Q 0.30 0.55 33° 25° 1.33/?* 4.60i?* 3 X 10- 8 M Q yr" 1 10 3 km s" 1 
model II (Casares et al. 2005) 8.2 R Q 0.23 0.72 70° 15° 1.32-R* 8.13.R* 10- 8 M Q yr" 1 10 3 km s" 1 



We perform the calculations of the synchrotron X-ray 
and the cascade IC 7-ray spectra for the location of the mag- 
netar at the periastron and the apoastron passages for both 
sets of parameters. The surface magnetic field of the massive 
star is investigated in the range B* = 10 2-4 G and the sur- 
face magnetic field of the magnetar is fixed to Bns = 3 x 10 14 
G. Electrons are accelerated in the transition region with 
the differential power law spectrum extending to the maxi- 
mum energies given by Eq. 1141 For the illustration we chose 
the spectral index equal to -2 and the acceleration parame- 
ter £_i = 1. The power in relativistic electrons is estimated 
from Eq. 1151 applying specific energy conversion efficiency r\. 
Note, that this power scales with the product of r\P^ 2 which 
is taken in our calculations to be equal to 10 -2 . We apply 
the inclination angles of the binary system to be consistent 
with other basic parameters (see Table 1) and assume the 
mass loss rate of the massive star equal to 3 x 1O _8 M0 yr -1 
(model I) and 10~ 8 Af© yr" 1 (model II). The results of the 
example calculations with the parameters mentioned above 
are shown in Fig. la-f. For all considered model parameters, 
the 7-ray spectra change drastically between the periastron 
and apoastron passages. No any emission above ~1 TeV is 
observed at the periastron due to the inefficient accelera- 
tion of electrons. The TeV 7-ray emission clearly extends 
through the TeV energy range at the apoastron, in agree- 
ment with the MAGIC and VERITAS observations. Note 
that 7-ray spectra depend also on the surface magnetic field 
of the massive star. For larger B*, secondary e pairs from 
the IC cascade lose more energy on the synchrotron process. 
That's way, the synchrotron X-ray spectra are on the higher 
level for larger B*. 

The distances of the periastron passage, r p , for both 
models are very similar (Table 1), but distances of the apoas- 
tron passage, r a , differ significantly. This can explain the 
stronger absorption features in the case of the IC cascade 
spectrum at the apoastron in the case of parameters pro- 
posed by Grundstrom et al. (2007) . However the steeper cas- 
cade IC spectra obtained at the periastron for the Casares et 
al. (2005) parameters has to be related to the differences in 
the angle of the periastron passage in respect to the location 
of the observer. Note, that in most of the considered cases 
(e.g. figures a,c,e,f), the GeV 7-ray fluxes does not change 
drastically between the periastron and apoastron passages. 
Therefore, it is likely that observed strong modulation of 
the TeV 7-ray emission with the period of the binary sys- 
tem may not be observed by the Fermi LAT detector in the 
GeV energies. 

The synchrotron spectra have more complicated shape 
due two the components produced by primary electrons and 
secondary e pairs. The strength of the component pro- 
duced by secondary e pairs strongly depends on the sur- 
face magnetic field of the massive star, B*, and on the dis- 
tribution of the magnetic field around the star. In general, 



for stronger B*, the synchrotron spectrum from secondary 
e pairs is stronger (see dashed curves in Fig. [1]), reaching 
saturation for B* ~ 10 3 G. At the apoastron passage, the 
second synchrotron bump becomes comparable to the first 
bump only for the largest considered values of the surface 
magnetic field of the massive star, i.e. B* = 10 4 G. 

Note that a part of primary electrons can likely escape 
from the transition region with the matter which can be 
expelled along the rotational axis of the magnetar, i.e. within 
the possible transient jet which is sometimes expected in 
the case of the accretion process occurring in the propeller 
regime. The synchrotron emission of these electrons in the 
jet is not taken into account in this calculations. 

In this simple model, the conversion coefficient of en- 
ergy from the transition region to relativistic electrons, 77, is 
kept independent on the phase of the binary system. How- 
ever, r\ may also depend on the accretion rate (i.e. a distance 
between the stars) in the way which is at present unknown. 
Therefore, any absolute normalization of the spectra is prob- 
lematic since it may vary with the phase of the binary sys- 
tem, e.g. as observed in the hard X-rays (Chernyakova et 
al. 2006). Only detailed phase dependent multi- wavelength 
observations can put constraints on the injection details of 
electrons at different phases of the binary system. 



5 CONCLUSION 

We show that strongly magnetized, accreting neutron stars 
(magnetars) close to the massive stars inside the binary sys- 
tems can be responsible for the acceleration of electrons up 
to TeV energies. These electrons interact with the magnetic 
field of the inner magnetosphere of the magnetar and with 
the soft radiation field from the massive companion produc- 
ing X-rays (synchrotron process) and 7-rays (inverse Comp- 
ton process). A part of injected 7-rays can be absorbed in the 
radiation of massive star developing IC e pair cascade. Sec- 
ondary e pairs lose significant part of their energy on the 
synchrotron emission. The 7-ray spectra produced in such 
scenario can extend up to TeV energies and the synchrotron 
spectra can extend up to MeV energies. The processes dis- 
cussed above concern to a specific stage in the evolution of 
the binary system, when the accretion rate onto the neutron 
star is at the range which guarantee the propeller phase. In 
fact, this phase may be short lived due to the variable wind 
from the massive star. When the accretion rate is different, 
the magnetar can move to the ejector or the accretor phase 
during which different scenarios should be considered. For 
example in the accretor phase the magnetar can gain an- 
gular momentum with the accreting matter and decrease 
its rotational period as expected in the case of the millisec- 
ond pulsars. We expect (based on the comparison of Eq. [16] 
and I17p that the propeller phase of the accretion onto the 
neutron star, in which observable fluxes of 7-rays can be 
detected, is relatively short lived. This may explain a small 
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Figure 1. The X-ray spectra (curves) and TeV 7-ray spectra (histograms) calculated for the periastron (thick) and apoastron (thin) 
passages of the magnetar inside binary system LSI +61 303 for two different sets of the parameters (see Table 1): figures (a), (b), and 
(c) for the model I and (d), (e), and (f) for the model II. The surface magnetic field of the massive star is = 10 2 G (a) and (d); 
B* = 10 3 G (b) and (e), 10 4 G (c) and (f). The calculations have been done for rip^ 2 = 10~ 2 and £_i = 1. For comparison we show 
the level of the X-ray and TeV 7-ray emission reported by RXTE (Smith ct al. 2008) and MAGIC (Albert et al. 2006) and VERITAS 
(Acciari et al. 2008) telescopes. 



number of massive binaries discovered up to now in the TeV 
7-rays. Applying the parameters of the binary system LSI 
+61 303, we show that such a model can explain the TeV 
7-ray emission observed close to the apoastron passage of 
the neutron star. The lack of observable TeV emission from 
the parts of the orbit close to the periastron is due to the 
inefficient acceleration of electrons and the suppression of 
the TeV 7-rays in the IC e cascade process in which also 
the synchrotron energy losses play an important role. 
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